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I discuss possible failures of local quark-hadron duality in the system of heavy-light quark mesons. As an 
example, I consider a correlator of two currents comprising heavy quark operators, and I compare the OPE 
expression with the result obtained by a complete insertion of hadronic states. After a smearing procedure, OPE 
and hadronic spectral functions agree with each other. However, the local behaviour of the two functions is 
different, and the difference manifests itself in a term which is absent in the OPE. 



1. Introduction 

Quark-gluon/hadron duality represents a ba- 
sic concept in the theoretical description of inclu- 
sive hadronic processes. In general, it means that 
some hadronic rates can be computed as the cor- 
responding partonic rates, and that high energy 
processes can be computed in terms of hadronic 
matrix elements of operators in an expansion, the 
Operator Product Expansion (OPE), whose lead- 
ing term is represented by the perturbative QCD 
expression. It is duality which allows us to ex- 
trapolate from the deep Euclidean region, where 
the OPE is defined, to the Minkowski domain, 
where physical observables are measured. 

An OPE-based approach has been recently pro- 
posed for the calculation of semileptonic and non- 
leptonic decay rates of hadrons containing one 
heavy quark M. The basic idea consists in the 
expansion in inverse powers of the heavy quark 
mass m,Q. In this approach, however, a distinc- 
tion must be maintained between semileptonic 
and nonleptonic decays. In general, hadronic 
and OPE amplitudes cannot be identical even at 
very high momentum transfer, due to the differ- 
ent production thresholds of multiparticles and 
of quarks and gluons. In the case of semilep- 
tonic heavy hadron decays this difficulty can be 
avoided, since one has to integrate over lepton 
variables, which amounts to a smearing of the 
OPE width. The equality between smeared OPE 
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and hadronic widths is referred to as global du- 
ality (different from local duality, i.e. without 
smearing). It is generally believed that global 
duality holds between quark-gluon and hadronic 
cross sections 0; for B and A^ semileptonic de- 
cays global duality has been extensively analyzed 
in §,§. 

For nonleptonic heavy hadron decays there are 
no lepton momenta to be integrated; therefore, 
one cannot invoke global duality to prove the 
identification of OPE and hadronic observables, 
and therefore local duality must be assumed. In 
this case, the validity of the OPE-based 1/ttiq 
expansion for the computation of nonleptonic 
heavy hadron inclusive decay rates appears to 
be more debatable. Indeed, in [0J it was sug- 
gested that the discrepancy between the predic- 
tion —7 — ^ > 0.9 101 and the experimental result 
r[B d ) " 

-7^4 = °- 78 + °- 04 might be solved assuming 

a violation of local quark-hadron duality, with the 
appearance of a 0(-^—) correction not predicted 
by OPE. 

At the moment, understanding the origin of 
this possible correction to the OPE is a difficult 
task. However, it is interesting to investigate pos- 
sible violations of local duality in some definite 
model. 

Several recent studies can be found in the 
literature concerning the validity of the quark- 
gluon/hadron duality in connection with non- 
perturbative QCD applications in the heavy 
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quark sector g-|12|. In particular, in Q the two- 
point function involving the difference between 
scalar and pseudoscalar heavy-light currents was 
considered as a benchmark to investigate global 
aa well as local duality properties. In the chiral 
limit (for the light flavour) and the infinite mass 
limit (for the heavy flavour), the coefficients of the 
OPE for the correlator can be calculated analyt- 
ically. Then, the spectral function derived form 
a particular hadronic model (in the time-like re- 
gion) is expanded in a power series (in the space- 
like region) and compared with the exact result. 
From this comparison one can gauge the validity 
of (local) duality. 

Here, I want to discuss how the results are sta- 
ble against modifications of the hadronic model. 
Evaluating the two-point function by a complete 
insertion of states derived from a relativistic con- 
stituent quark model, I'll show that an explicit 
violation of local duality can be detected, in the 
form of an unexpected correction to the OPE JToJ . 

2. Example of violation of local duality 

I consider the correlator (8j 

11(g) = i / dxe«*[< 0\T(J s (x)4(0))\0 > - 



< 0|T(J P (x)4(0))|0 >] = U s (q) - U P (q){l) 

with J s (x) = Q(x)q(x), Jp{x) = Q(x)i-f 5 q(x), 
Q(x) and q(x) being heavy and light quark oper- 
ators, respectively. In the chiral limit, m q — > 0, 
IT(<7) vanishes in perturbation theory. In the in- 
finite heavy quark mass limit, toq — > oo, it is 
convenient to choose q^ — (nig — e, 0), so that 
the correlator becomes a function of e: 

11(e) = l£°°dre-^ $(r) (e > 0) , (2) 

$(t) representing the non-local quark condensate 

$(r) = <0\q(T)e ig °f° dy " A » iy) q(Q)\Q> . (3) 
When e >> Aqcd, the OPE expression for 11(e) 

< o\qq\o > „ 



n, 



OPE 
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9 4 fi 
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is given in terms of the quark condensate 
< 0\qq\0 >= (-240 MeVf and of the 
mixed quark-gluon condensate, parametrized by 

ml = < > = .8±0.2GeV». 

< 0\qq\0 > 

The positive coefficients C2„ depend on the actual 
form of the non-local condensate; notice explicitly 
the alternating signs in Eq.(Q) and the absence of 
even powers of e _1 . 

It is possible to identify a model for 11(e) having 
the expansion (||) [||: 
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A 



4A 



2A 



with A a parameter and 



z + l 
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(5) 
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ip(z) being the logarithmic derivative of the 
Gamma function. 11(e) in (m can be written as 



11(e) 



< Q|gg|Q > 

2A 



E 



c-i)j 



U e/A + 2j + l 



and admits the asymptotic expansion 
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< Q|gg|Q > 
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(Em Euler numbers). Comparison of (||) with (^) 
indicates that the model is able to reproduce the 
right power structure of 11(e). The spectral den- 
sity associated to (|j) displays an infinite number 
of equally spaced poles located along the negative 
e axis, with residues having alternating signs. 

The correspondence between (||) and (^), how- 
ever, follows from the assumed form (|5|) for the 
correlator. It would be interesting to investigate 
the properties of global and local duality using a 
representation of the correlator obtained from a 
hadronic model more closely related to QCD. 

Let us compute ([!]) by a complete insertion 
of hadronic states. We consider the N c — > oo 
limit, since in this case the only contributions 
are J p = + and CP single particle states, con- 
tributing respectively to lis and lip. We denote 
by \S n > and \P n > the states, with masses 
Ms n and Mp n respectively, and we define the 
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current-particle matrix elements: < OjJslSVi >= 



Mr. 



Mi 



m ffs n , < 0\J P \P n >= ^/p„. In the limit 
m,Q — > oo the masses can be written as: 



M 



S„(P„) 



0(— ) 



(9) 



The binding energies <5s„(p„) can be obtained by 
solving the relativistic wave equation 



-V 2 



-V 2 



■ ml + V(r)]V n (r) = M„*»(rfl0) 



where the (central) potential V(r) is chosen as 
V(r) = fj, 2 r, with constant [i (string tension). 
Masses, wavefunctions and current-particle cou- 
plings of scalar and pseudoscalar states can be 
computed applying WKB methods to ( |l(i| ) jL3l 
Pf. One gets: 



An) 



An) 



ln{2n + l+-) 



(n) 



(11) 



(12) 



Assuming the WKB approximation for the 
lowest-lying state, an estimate can be derived for 
the string tension: for mj, = 4.6 — 4.7 GeV one 
has /i ~ 300 MeV. 

The choice of considering only single particle 
states between the currents in the correlator is 
not too restrictive, as it can be shown by con- 
sidering the imaginary part of the correlator of 
scalar currents lis (lip) in Eq.((l|). Computing 
the imaginary part of the quark loop diagram one 
has, for E = — e >^ +oo: 

(E) - ^ ■ (13) 



ImW° s PE 



8tt 



On the other hand, the resonance model gives: 

had _ 



ImIL 



71=0 



§\E- n^2w(n + 7/4) .(14) 



Although the two expressions (|l|) and ( pi) look 
very different, after a smearing procedure 2] ob- 
tained by J2 n ~ * / one nas 



Im n 



had 
S 



dn5 £'-Aiv / 27r(n + 7/4) 



3E 2 
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= ImU 



OPE 



(E) 



(15) 



Therefore lis(E) and Up(E) satisfy global dual- 
ity, at least at the leading order for E — > oo. 

Let us now compute the correlator (Q) by the 
hadronic insertion. We get 



q 3 00 
U had (e) = ^ ^ 



(-l)"+Vn + 3/2 
8 ^ n^o e + MvM™ + 3 / 2 ) 
and (E > 0) 

R had {E) = ImU had (E) 
3fi 2 E 



(16) 
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.(17) 



From the first four terms in (|4[) we have: 
R° PE (E) = Im U OPE (E) 

2 

= -I < 0\q q \0 > [S(E)-^5"{E) 

+ c ^v) { E)-c^ VI \E)..] , (18) 

• ! , 5 61, _ 

with unknown C4, cr (04 = — , Cr = m the 

4, b v 4 64 , 

model H). Also in this case the two expressions 
look very different, but they can be compared af- 
ter an appropriate smearing; following (2) we con- 
sider the smeared function 



A) 
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dE' 



R had (E') 



(E — E' 



A 2 



(19) 



and a similar expression for R OPE (E, A). For 
A —> on can prove that R(E, A) -» i?(£). 
Moreover, since 

R(E, A) = — [Il(E - iA) - U(E + iA)] , (20) 

for A >> Aqcd one is far from physical singu- 
larities. 

For a large energy E the two expressions R OPE 
and j?' md should coincide regardless of A. Such 
a request implies: 



< 0|gg|0 >= 



0.51 /x d 



(21) 



e.g. n = 317 MeV, in agreement with the value 
derived from the spectrum. For smaller values of 
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E, i.e. 1 < E < 20 GeV, the numerical results 
concerning R had (E, A) and R OPE (E 1 A) are re- 
ported in Fig.l, where the ratio P = R had /R OPE 
is plotted as a function of E for several values of 
A. As expected, the agreement between R had and 
R OPE improves with increasing A; in particular, 
for A = 3.0 GeV the difference does not exceed 
20%. For any value of the smearing parameter 
A the ratio P tends to unity for large energies. 
Therefore, the resonance model satisfies the re- 




Figure 1. Plot of the ratio P = R had /R OPE . 
From top to bottom: A = 1.5, 2.0, 3.0, 5.0 GeV 



quirement of global duality: the smeared imagi- 
nary parts of the correlator, when computed by 
OPE or by hadronic states, agree with each other. 

Concerning local duality, the expressions are 
significantly different, since 



rr ad ( e ) 



0.5V 1 



1 



mo 
e 



0.93- 



.(22) 



The factor multiplying e 1 coincides numerically 



with 



<0|gg|0> 



for (i = 317 MeV (eq.(Elh). As 



for the mass parameter mo, numerically we find 
mo = 560 MeV. In conclusion one has: 



U had (e) 



< 0|gg|0 > 
4e 



1 



m 
e 



™2 



0.93- 



.(23) 



The comparison between (£3|) and (Q) shows a 
violation of local duality, which manifests itself in 
the form of an unexpected term in the asymptotic 

expansion in powers of — . 



3. Conclusions 

In the case of a simple correlator of quark cur- 
rents, using a particular quark model and in the 
N c — > oo limit, I have shown an example of viola- 
tion of local duality which occurs in a condition 
where global duality is verified. The difference 
between the hadronic and the OPE spectral func- 
tions is made evident by a term which is absent 
in the expansion predicted by OPE. 

Of course, the calculation of the correlators re- 
quired for the evaluation of the Bd and Ah inclu- 
sive decay widths, is, beyond our present possi- 
bilities. However, the example supports the con- 
jecture that the A;, lifetime problem can be ex- 
plained by the presence of a 1 /tuq correction not 
included in the usual OPE expansion. 
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my opinion, is instructive is to see how global du- 
ality is recovered, in terms of the energy E and of 
the smearing parameter A. 



Discussions 



B. Blok (Technion, Haifa) 

How are your results model dependent? 

P. Colangelo 

As discussed, there is a dramatic difference be- 
tween two choices of the hadronic representations 
of the correlator Therefore, the results are 

model dependent, although the model I presented 
seems rather realistic and QCD-inspired. The 
main point is that it allows to detect a mechanism 
of violation of local quark-hadron duality which is 
similar to what is needed to explain the A^ life- 
time problem. 



M. Neubert (CERN, Geneva) 
In my opinion, the violation of global duality 
found in your model (~ 10— 40%) are surprisingly 
large. Is the E-dependence of these violations in 
accordance of general expectations as discussed, 
e.g., in the previous talk by B.Blok? 

P. Colangelo 

The numerical details can be ascribed to some fea- 
tures of the model, for example to the assumption 
of vanishing width of the resonances. What, in 



